Cupric chloride C11CI2 as an S = 1/2 chain multiferroic 
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Magnetoelectric properties were investigated for an S = 1/2 chain antiferromagnet CuCb, which 
turns out to be the first example of non-chalcogen based spiral-spin induced multiferroics. Upon the 
onset of helimagnetic order propagating along the fr-axis under zero magnetic field (if), we found 
emergence of ferroelectric polarization along the c-axis. Application of H along the 6-axis leads to 
spin-flop transition coupled with drastic suppression of ferroelectricity, and rotation of H around the 
6-axis induces the rotation of spin-spiral plane and associated polarization direction. These behaviors 
are explained well within the framework of the inverse Dzyaloshinskii-Moriya model, suggesting the 
robustness of this magnetoelectric coupling mechanism even under the strong quantum fluctuation. 

PACS numbers: 75.85.+t, 77.80.Fm, 75.45.+j 



I. INTRODUCTION 

Multiferroics, materials with both magnetic and dielec- 
tric orders, have attracted revived interest. 1 While the 
coupling between these orders are weak in general, re- 
cent discovery of spin-driven ferroelectricity in frustrated 
helimagnets has enabled unprecedentedly large magneto- 
electric (ME) effects such as flopj 2 -"— reversal^ or rotation^ 
of electric polarization (P) under applied magnetic field 
(H). Here, the key problem is the coupling mechanism 
between ferroelectricity and helimagnetism. The most 
successful scheme to explain such ME coupling is the 
inverse Dzyaloshinskii-Moriya (IDM) model, 7 which de- 
scribes the local polarization pij produced between two 
magnetic moments Si and Sj as: 



Pij 



Aeij x (Si x Sj)^ 



(i) 



where e^- is an unit vector connecting two magnetic sites, 
and A is a coupling coefficient related to the spin-orbit 
interaction. Since the vector spin chirality (Si x Sj) is 
perpendicular to the spin-spiral plane, this model pre- 
dicts that the ii"-induced tilt of spin-spiral plane leads 
to directional change of P. Ferroelectric (FE) and ME 
natures in several classical helimagnets like RMnOsr 1 ^ 
NisV20^ and MnW04 4 are successfully explained by 
this IDM scheme. 

In contrast, still in controversy is the ME coupling 
mechanism in S = 1/2 chain magnets, where strong 
quantum fluctuation is believed to have some profound 
effects on their ME response^^ Typical examples are 
LiCu202 and LiCuV04, both of which are characterized 
by edge-shared Cu02 chains. While their simple crys- 
tal structures and reported helimagnetism are seemingly 
typical of the IDM scheme, the experimentally observed 
P-direction or ME response for applied H appears to 
contradict with its prediction. For example, LiCu2 02 



hosts P || c in the helimagnetic ground stated but its 
proposed magnetic structures are contradictory among 
several experiments i^r— Even if we simply assume bc- 
spin spiral consistent with the IDM model, P-behaviors 
under applied H contradict with its naive prediction^ 
In case of LiCuV04, P || a is found in afr-spiral spin 
phase at T, consistent with the IDM mode h 16 i 17 How- 
ever, the observation of P || c, assigned to the frospiral 
spin phase, 18 was not reproduced by another group 
Since both compounds frequently contains Li-Cu inter- 
substitution due to their close ionic radii, Moskvin et 
al. have claimed that the observed FE and ME natures 
stem purely from exchange striction and crystallographic 
defects, not from the spin-orbit interaction (or the IDM 
scheme) Furthermore, some recent theoretical study 
predicted that quantum fluctuation may largely reduce 
the effective magnitude of P induced via the spin-orbit 
interaction^ To testify the validity of the IDM model 
in quantum chain magnets, it is crucial to check the ME 
response in other 5=1/2 compounds with similar edge- 
shared chain structures. 

Our target compound, anhydrous cupric chloride 
CuCl2 crystalizes into distorted Cdl2 form with mono- 
clinic C2/m space group and j3 = 122°. 23 While origi- 
nal Cdl2 structure consists of the stacking of triangular 
lattices along the z-axisj2i they are extended along the 
a- axis due to Jahn- Teller active Cu 2+ ions (Fig. DJa)). 
As a result, CuCl2 can be regarded as the aggregate of 
edge-shared chains running along the 6-axis, with CuCU 
square plaquettes lying in the 6c-plane (Fig. []Jb)). Mag- 
netism is dominated by the intra-chain coupling between 
Cu 2+ (S = 1/2) ions, and competition between ferromag- 
netic nearest-neighbor interaction and antiferromagnetic 
next-nearest-neighbor interaction stabilizes the helimag- 
netic ground state below 24 K^~— Recent powder neu- 
tron scattering study suggested the cycloidal magnetic 
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FIG. 1: (color online), (a)-(c) Crystal structure of CuCb, and 
P-direction observed at the ground state. The fcc-cycloidal 
spin order suggested by Banks et a/.— is illustrated in (b), 
and also in (c) with solid rounded square representing spin- 
spiral plane. Dashed rounded square indicates the possible 
tilting of spin-spiral plane as revealed in this study (see text), 
(d) ESR signal taken at room temperature under various di- 
rections of if confined within the ac-plane. Oh is defined as 
the angle between the a-axis and if-direction. Each dashed 
line represents a fitted curve with a single Lorentzian reso- 
nance, (e) Angle dependence of ^-factor. 



order propagating along the b- axis, with spin spiral con- 
fined in the frc-plane (Fig. DJb)) and propagation vector 
q ~ (1, 0.226, 0.5). 27 While no dielectric measurements 
have been reported, the latest calculation based on den- 
sity functional theory (DFT) predicts emergence of ferro- 
electricity along the c-axis. 27 In this paper, we report the 
experimental discovery of FE and ME natures in CuCl2, 
and prove that the IDM mechanism is still robust even 
under the strong quantum fluctuation. CuCl2 is also the 
first example of non-chalcogen based spiral-spin induced 
multiferroics, which promises further discovery of unique 
ME function in many MJ\~2-type compounds and other 
forms of halide compounds. 



II. EXPERIMENT 

Single crystals of CuCl2 were grown by a Bridgman 
method. They were cleaved along planes perpendicular 
to the z-axis, and cut into a rectangular shape with addi- 
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FIG. 2: (color online). Temperature dependence of (a) mag- 
netic susceptibility (b) dielectric constant e, and (c) electric 
polarization P. In (d) and (e), if -dependences of magnetiza- 
tion M, e, and P under if || b are indicated. Large and small 
arrows denote corresponding ordinate scale of physical quan- 
tity and direction of field scans, respectively, (f) H-T phase 
diagram for if || b. Circles, squares, triangles are the data 
points obtained from M, e and P profiles, respectively. Open 
(closed) symbols are taken from T- (if-) increasing runs. 



tional faces perpendicular to the a- or 6-axis. Silver paste 
was painted on end surfaces as electrodes. Due to its 
moisture sensitivity, the specimen was handled in an Ar- 
filled glove box. To deduce P, we measured the polariza- 
tion current with constant rates of temperature(T)-sweep 
(5 ~ 20 K/min), if -sweep (100 Oe/sec), or if -rotation 
(2° /sec), and integrated it with time. To enlarge the 
population of specific P-domain, the poling electric field 
(E = 150 ~ 400 kV/m) was applied in the cooling pro- 
cess and removed just prior to the measurements of polar- 
ization current. Dielectric constant e was measured at 1 
MHz using LCR- meter. Magnetization M was measured 
with a SQUID magnetometer. ESR signal was measured 
by JEOL JES-FA200 at X-band frequency (~ 9.0 GHz). 



III. RESULTS AND DISCUSSIONS 

As suggested in Ref. 27, the ac-twin domains are ex- 
pected to readily occur. To check this possibility, we first 
performed ESR measurements under various directions 
of if confined within the ac-plane(Fig. [T](d)). Hereafter, 
we define Oh as the angle between the a-axis and in- 
direction. Each observed profile can be fitted well with 
a single Lorentzian resonance for all Oh, indicating our 
crystal grown by the Bridgman method has no crystallo- 
graphic twinning. The deduced ^-factor shows sinusoidal 
^-dependence (Fig. [He)), whose maximum and mini- 
mum values agree well with those previously reported 
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FIG. 3: (color online). P-E hysteresis loop for CuCl2 mea- 
sured with E || z. Electric field was swept at the rate of 
13kV/m-sec. 



for a twinned crystal^ 

Next, we measured T-dependence of magnetic suscep- 
tibility Xi e ? an d P (Figs. [2ja)-(c)). x suddenly drops 
at Tn ~ 24 K, which signals the transition into a spi- 
ral magnetic phase. Simultaneously, z-component of e 
(e z ) shows a sharp anomaly and a- and ^-components 
of P (P a and P z ) begin to develop. The direction of 
P z can be reversed with reversal of applied E, produc- 
ing the typical P-E hysteresis curve (Fig. [3]), while no 
P\) component could be confirmed. These results imply 
strong correlation between helimagnetic and FE orders 
in CuCl2. Based on the frc-plane helimagnetic structure 
suggested in Ref. 27, the IDM model as well as the DFT 
calculation^ predicts \P a /P z \ ~ 0.64 (i.e. P \\ c). This 
roughly agrees with the observed \P a /P z \ ~ 0.70. 

Figures [2fd) and (e) indicate independence of M, e, 
and P for H || b. At 4 T, M-profile shows a clear step as 
already reported ] 27 ! 28 Concomitantly, e z shows a sharp 
peak and both P a and P z are drastically suppressed. 
Since antiferromagnetic spins favor to lie within a plane 
perpendicular to H, this transition should correspond to 
a spin-flop into the ac-spiral spin state. The ac-spiral spin 
structure belongs to a magnetic form called proper screw, 
where spin-spiral plane is perpendicular to the modula- 
tion vector along the 6-axis. The IDM model predicts 
P = for this type of spin order due to the relationship 
(Si x Sj) || eij || 6, which is consistent with the observed 
suppression of P. Figure [2ff) summarizes the obtained 
H-T phase diagram for H || b. The boundary of the FE 
phase always coincides with that for magnetic phases, 
which proves the interplay between FE and helimagnetic 
natures. 

We further investigated the properties under H JL b. 
Here, we adopt the same definition of Oh as used for ESR 
measurements. Figure[H(a) indicates iJ-dependence of M 
measured at various Oh- While no magnetic transition 
has been reported for H JL 6, 28 we found a clear signa- 
ture of spin-flop at H$f ^ 4 T most pronounced around 
Oh = 100°. ^-dependence of x(= M/H) was also mea- 
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FIG. 4: (color online), (a) if -dependence of M taken under 
various directions of H around the 6-axis. The base lines 
of data are arbitrarily shifted, (b) Angle-dependence of x 
measured under H JL b. Dashed lines indicate the fits with 
sinusoidal function, (c)-(e) Corresponding development of z- 
component of e (e z ) as well as z- and a-components of P (P z 
and P a ). Before measurements of P, the specimen was cooled 
at Oh = with poling E applied along the z-axis. Dashed 
lines indicate the behavior expected from Eq. (2). Arrows 
indicate the direction of if-rotation. 



sured (Fig. [H(b)), and x sinusoidally changes with mini- 
mum at Oh ~ 100° below Hsf- In general, the sharpest 
transition of spin-flop as well as the minimum value of 
X should be observed when H is applied parallel to the 
magnetic easy-plane. These results imply the magnetic 
easy-plane, i.e. spin-spiral plane at the ground state, is 
tilted from the originally suggested frc-plane toward the 
6z-plane by about 20° (Fig. [He)). Above H SF , x still 
modulates sinusoidally but with different x-minimum po- 
sition at H ~ 122° (i.e. H \\ c). With H > H S f, the gain 
of Zeeman energy exceeds the energy loss due to magnetic 
anisotropy, and continuous rotation of spin-spiral plane 
is expected. In this case, ^-dependence of x rather re- 
flects the anisotropy of g- value, 29 whose minimum is also 
confirmed to appear at H \\ c (Fig. [2(e)). 

To investigate the behavior of P under H rotating 
around the 6-axis, we simultaneously measured P z and 
P a using two pairs of electrodes. Thus, both P and H 
can be expressed as vectors within the ac-plane. We also 
define as the angle between the a- axis and observed 
P-direction (Fig. [5](d)). Figures IH(d) and (e) indicate 
^-dependences of P z and P a measured at 5 T. When H 
is rotated by 180°, P-direction is always found to be re- 
versed. To see the behavior of P more straightforwardly, 
the trace of P is plotted in the P a -P z plane (Fig. [51(a)). 
It forms a shape like elongated ellipse. In Figs. 0(b) 
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FIG. 5: (color online), (a) Trace of P a and P 2 under if 
rotating around the 6-axis. (b),(c) Magnitude and direction 
of P as a function of if-angle. Arrows indicate the direction 
of i7-rotation. The data are taken from Figs. 5] (d) and (e), 
and the setup for measurements is shown in (d). Op (Oh) 
is defined as the angle between P- (H-) direction and the 
a-axis. (e),(f) The expected relationship between P, H and 
spin-spiral plane (depicted as rounded square). 



and (c), ^-dependences of \P\ (magnitude of P) and 
Op are indicated. Op takes almost constant value around 
Op = 120° or 300°, suggesting the major axis of observed 
P-ellipse is pointing at the c-axis. If we assume that 
H is always perpendicular to the spin-spiral plane, i.e. 
H || (Si x Sj), the IDM model predicts P = P 1 || c for 
H -L c (Fig. E(e)) and P = P 2 _L c for H \\ c (Fig. EJf)). 
For general Oh, P is given as 

P = Pi sin(122° - H ) + P 2 cos(122° - H ), (2) 

which forms an ellipse-shaped trace with Pi and P 2 as the 
major and minor axes, respectively. From the |P|-profile, 
we deduced |Pi| ~ 31 /iC/m 2 and |P 2 | ~ 2 /iC/m 2 . 

In Figs. Hfd) and (e), the P-behavior expected from 
Eq. (2) is plotted as dashed lines. While the calculated 
P-profile roughly agrees with the observed one, small gap 
still exists between them. This deviation reverses its sign 
at Oh ~ 100°, where H becomes parallel to the mag- 
netic easy-plane. Correspondingly, e also shows small 
anomaly at Oh ~ 100° (Fig. Sfc)). These behaviors can 
be well explained by assuming that the spin-spiral plane 
is tilted from the originally expected H \\ (Si x Sj) posi- 
tion toward the magnetic easy-plane. A similar effect of 
magnetic-anisotropy drag on P has also been observed in 
the iJ-rotating experiment on Eui-^Y^MnOa. 6 

Thus, we conclude that the IDM scheme can well 
reproduce the observed FE and ME natures, even for 



CuCl 2 with 5 = 1/2 quantum spin chains. Note that 
P || c relationship observed at T can be justified even 
with a slight revision of the originally suggested frospiral 
spin structure, since deduced ratio |Pi|/|P 2 | ~ 15 is 
quite large. Notably, when H and spin-spiral plane is 
rotated counter-clockwise, P is found to rotate clock- 
wise (Fig. |5(a)). This is in contrast with the case for 
Eui-^Y^MnOa, 6 where both H and P rotate in the same 
direction. The observed manner of P-rotation and large 
lAl/l^l ratio are in accord with the recent DFT calcula- 
tion for edge-shared Cu0 2 chain compounds, 30 and these 
features would reflect the anisotropy and sign of coupling 
coefficient A in Eq. (1). 

CuCl 2 is also the first example of non-chalcogen based 
spiral-spin induced multiferroics. While observed |Pi| ~ 
31 /iC/m 2 is somewhat smaller than the calculation 
(|Pi| ~ 84 /iC/m 2 ), 27 it is comparable with the case for 
other helimagnetic oxides (2000 ~ 5 /iC/m 2 ). Interest- 
ingly, a recent theoretical study suggested the choice of 
anion may largely affect the value of induced P through 
the different strength of metal-ligand hybridization and 
spin-orbit coupling. 31 This means that the ME response 
can be enhanced if we choose appropriate anion as the 
ligand. Until now, the study of ferroelectric helimagnets 
is almost limited to oxides, partly because their isostruc- 
tural chalcogen relatives with larger anions (i.e. sulfides 
or selenides) are often electrically too leaky to perform 
dielectric measurements^ 2 -"— Since halogens have larger 
electronegativity than chalcogens, halides are better in- 
sulating and enable the investigation of ME properties for 
a wider variety of anions. For example, most of MX 2 - 
type halides with X = CI, Br, and I consist of stacking 
of undistorted triangular lattices, which realizes various 
types of spiral spin orders^ - — The systematic investi- 
gation of FE properties of whole MX 2 system will of- 
fer a good opportunity to clarify the anion dependence 
of magnetically-induced ferroelectricity, which may con- 
tribute to a general strategy to obtain larger ME re- 
sponses. 



IV. CONCLUSIONS 



In summary, we have experimentally revealed 
magnetically-driven ferroelectricity in an S = 1/2 chain 
helimagnet CuCl 2 . Observed P-behaviors under applied 
H can be reproduced well within the framework of the 
inverse Dzyaloshinskii-Moriya model, suggesting the ro- 
bustness of this ME coupling mechanism even under the 
effect of strong quantum fluctuation. CuCl 2 is the first 
example of non-chalcogen based spiral-spin induced mul- 
tiferroics, which promises further discovery of unique 
magnetoelectric function in many MX 2 -type compounds 
and other forms of halide compounds. 
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